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Abstract: We study leptogenesis in the supersymmetric standard model plus the seesaw. 
We identify important qualitative differences that characterize supersymmetric leptoge- 
nesis with respect to the non-supersymmetric case. The lepton number asymmetries in 
fermions and scalars do not equilibrate, and are related via a non-vanishing gaugino chem- 
ical potential. Due to the presence of new anomalous symmetries, electroweak sphalerons 
couple to winos and higgsinos, and QCD sphalerons couple to gluinos, thus modifying the 
corresponding chemical equilibrium conditions. A new constraint on particles chemical 
potentials corresponding to an exactly conserved -R-charge, that also involves the number 
density asymmetry of the heavy sneutrinos, appears. These new ingredients determine the 
3x4 matrices that mix up the density asymmetries of the lepton flavours and of the heavy 
sneutrinos. We explain why in all temperature ranges the particle thermodynamic system 
is characterized by the same number of independent quantities. Numerical differences with 
respect to usual treatment remain at the 0{1) level. 
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1. Introduction 

Leptogenesis [1,2] is a theoretical mechanism that can explain the observed matter-antimatter 
asymmetry of the Universe. An initial lepton asymmetry is generated in the out-of- 
equilibrium decays of heavy singlet Majorana neutrinos, and is then partially converted 
to a baryon asymmetry by anomalous sphaleron interactions [3]. Heavy Majorana singlet 
neutrinos are also a fundamental ingredient of the seesaw model [4] that explains in an 
elegant way the suppression of the neutrino mass scale with respect to all other particle 
masses of the Standard Model (SM). 

The discovery of neutrino oscillations has promoted leptogenesis to an utmost at- 
tractive scenario to explain the origin of the cosmic baryon asymmetry. This is because 
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without any fine-tuning of tlie seesaw parameters, a neutrino mass scale naturally com- 
patible with the solar and atmospheric neutrino mass square differences would be optimal 
to yield the correct value of the baryon asymmetry. The possibility of explaining two ap- 
parently unrelated experimental facts (neutrino oscillations and the baryon asymmetry) 
within a single framework has boosted the interest in leptogenesis studies, leading to im- 
portant developments in the field, as for example the inclusion of thermal corrections [5], 
spectator processes [6,7], flavour effects [8-12], CP asymmetries in scatterings [13], lepton 
asymmetries from the decays of the heavier Majorana neutrinos [14, 15], and many more. 

In spite of all these advancements, we believe that a proper treatment of leptogenesis 
in the supersymmetric case is still lacking. Supersymmetric leptogenesis constitutes a 
theoretically appealing generalization of leptogenesis for the following reason: while the 
SM equipped with the seesaw provides the simplest way to realize leptogenesis, such a 
framework is plagued by an unpleasant fine-tuning problem. For a non degenerate spectrum 
of heavy Majorana neutrinos, successful leptogenesis requires generically a scale for the 
singlet neutrino masses that is much larger than the electroweak (EW) scale [16] but 
at the quantum level the gap between these two scales becomes unstable. Low-energy 
supersymmetry (SUSY) can naturally stabilize the required hierarchy, and this provides a 
sounded motivation for studying leptogenesis in the framework of the supersymmetrized 
version of the seesaw mechanism.* Supersymmetric leptogenesis has been studied in several 
places, both in dedicated studies [18] or in conjunction with SM leptogenesis [5]. However, 
several features that are specific of supersymmetry in the high temperature regime relevant 
for leptogenesis, in which soft supersymmetry breaking parameters can be effectively set 
to zero, have been overlooked or neglected in these studies. When the new ingredients 
are left out, in spite of the large amount of new reactions, the differences between SM 
and supersymmetric leptogenesis can be resumed by means of simple counting of a few 
numerical factors [2, 19,20], like for example the number of relativistic degrees of freedom 
in the thermal bath, the number of loop diagrams contributing to the CP asymmetries, 
the multiplicities of the final states in the decays of the heavy neutrinos and sneutrinos. 

In this paper we show that, in contrast to the naive picture, supersymmetric lepto- 
genesis is rich of new and non-trivial features, and genuinely different from the simpler 
realization within the SM. A first important effect that follows from the requirement that 
the supersymmetry breaking scale should not exceed by much the 1 TeV scale, is that above 
a temperature T ~ 5 x 10^ GeV the particle and superparticle leptonic density asymmetries 
do not equilibrate. It is then mandatory to account in the Boltzmann equations for the 
differences in the number density asymmetries of the boson and fermion degrees of freedom, 
that can be given in terms of a non-vanishing gaugino chemical potential. A second feature 
is that when soft supersymmetry breaking parameters are neglected, additional anomalous 
global symmetries that involve both SU{2) and SU{3) fermion representations emerge [21]. 
As a consequence, the EW and QCD sphaleron equilibrium conditions are modified with 
respect to the SM, and this yields a different pattern of sphaleron induced lepton-flavour 

*In turn, supersymmetric leptogenesis introduces a certain conflict between the gravitino bound on the 

reheat temperature and the thermal production of the heavy singlets neutrinos [17]. In this paper we will 
not be concerned with the gravitino problem, nor with its possible ways out. 
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mixing [8,10,11]. In addition, a new anomaly-free i?-symmetry can be defined and the cor- 
responding charge, being exactly conserved, provides a constraint on the particles density 
asymmetries that is not present in the SM. Finally, a careful counting of the number of 
constraining conditions versus the overall number of particle density asymmetries reveals 
that four independent quantities, rather than the three of the SM case, are required to give 
a complete description of the various particle asymmetries in the thermal bath, with the 
additional quantity corresponding to the number density asymmetry of the heavy scalar 
neutrinos. In spite of all these qualitative differences, numerical corrections with respect 
to the case when all new effects are neglected remain at the 0(1) level. This is because 
only spectator processes get affected, while the overall amount of CP violation driving 
leptogenesis remains the same than in previous treatments. 

In Section || we start with a general description of the consequences of having different 
reactions dropping out of equilibrium as we imagine to raise the thermal bath temperature. 
We will next discuss in Section |2.2| the superequilibration regime, in which leptogenesis 
would proceed much alike in the SM case, but that in fact can be realized only at tem- 
peratures much below the lower limit for successful leptogenesis. In Section we discuss 
the most relevant theoretical issues, that are particle-sparticle non-superequilibration, the 
new global symmetries, the modified sphaleron conditions, and the relevant conservation 
laws that constrain the particle density asymmetries. We also analyze specific temperature 
ranges in which leptogenesis can be successful, and for each range we present results for 
the matrices that control the lepton flavour mixing induced by the fast sphaleron processes 
and that, differently from the SM case, now have four columns. In Section ^ we discuss the 
numerical relevance of our study with respect to previous treatments. Finally, in Section ^ 
we recap our main results and draw the conclusions. 

2. Chemical equilibrium conditions and conservation laws 

At each specific temperature, particle reactions in the early Universe can be classified 
according to their (thermally averaged) rates as: 

I Much faster than the Universe expansion rate. In the whole temperature range that 
we will consider, the top-quark Yukawa interactions and the SU{3) x SU{2) x U{1) 
gauge interactions are examples of reactions of this type. 

II Much slower than the expansion. For example, at temperature T ^ (1+tan^ /3)10^ GeV 
the rates of processes that involve the electron Yukawa coupling are completely neg- 
ligible with respect to the Universe expansion rate. 

Ill Of the order of the expansion. This class includes the neutrino Yukawa interac- 
tions responsible for leptogenesis when the temperature is of the order of the heavy 
Majorana neutrino mass M. 

Reactions of type I enforce specific conditions on the chemical potentials of the particles 
in the thermal bath. Reactions of type II imply that in the relevant temperature interval. 
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some Lagrangian parameters can be effectively set to zero. In most cases (but not in all 
cases) tliis results in exact global symmetries that correspond to conservation laws for the 
corresponding charges. These conservation laws constrain the asymmetries in the number 
densities of the particles that carry the conserved charges. For particles that are in chemical 
equilibrium, there is a statistical correspondence between the boson (6) and fermion (/) 
number density asymmetries An;, j = n^j — fibj and the corresponding chemical potentials 
//bj. In the relativistic limit rribj ^ T and at first order in /i^j <^ T this correspondence 
has the particularly simple form: 

Anfe = |rV6, Anj = ^T^f,f , (2.1) 

where Qhj are the number of degrees of freedom of the corresponding particles. Therefore, 
if chemical equilibrium is enforced for all particles entering a conservation law condition, 
the constraints on number density asymmetries implied by reactions of type II can be 
directly translated into constraints for the particles chemical potentials. However, when 
a charge conservation condition involves particle species that in general remain out of 
chemical equilibrium, as is the case for states that participate only in reactions of type III, 
such a correspondence does not exist, and the constraints must be formulated in terms of 
number density asymmetries. 

Conservation conditions related to gauge symmetries, like hypercharge conservation 



(see Eqs.(2/7)-(2^) below), are somewhat different from the conservation laws stemming 
from reactions of type II. This is because EW symmetry restoration at high temperature 
occurs when the order parameter v is dynamically driven to zero, which is different from 
the cases when some Lagrangian parameters can be neglected. However, for our purposes 
this difference has no relevance, and we will treat hypercharge conservation in the same 
way than other global conservation laws. 

Finally, reactions of type HI violate some symmetries, and thus spoil the corresponding 
conservation conditions, without being fast enough to enforce chemical equilibrium. To this 
class belong, by assumption, all the reactions that involve the neutrinos N, the sneutrinos 
N and the anti-sneutrinos A^*. The evolution of their number densities and, eventually, 
of the corresponding non-conserved charges, must then be tracked by means of Boltzmann 
equations. 

In the following we will analyze the set of conditions that constrain the particle number 
density asymmetries. Whenever possible we will express these constraints in terms of 
particles chemical potentials. This is the appropriate language for the fast reactions of 
type I, and in most cases can be consistently used also for conditions related to reactions 
of type II. In principle there are as many chemical potentials, or more precisely as many 
density asymmetries as there are particles in the thermal bath. However, there is also a 
large set of conditions enforced by processes in chemical equilibrium and by conservation 
laws which ensure that just a few independent quantities are sufficient to determine all the 
others. In carrying out our analysis we imagine to start from relatively low temperatures, 
that is, from temperatures well below the lower limit for successful leptogenesis. As we 
raise the temperature, both the expansion rate and the particle reactions rates vary, and 
this changes the way particle reactions are assigned to the previous three classes. More 
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specifically, we will only consider what happens when some B — L conserving reaction moves 
instantaneously from class I to class II, since for these reactions the regime in which the rates 
become of the order of the expansion (class III) is in general a short-lasting transient with 
no relevant effects associated to it. Since we are interested in leptogenesis, we also assume 
that, for each temperature regime, the heavy Majorana neutrino masses and couplings 
are such that B — L violating processes always belong to the third class. ^ According to 
the classification scheme described above, we would expect that at each temperature the 
thermodynamic system of the particle soup can be described in terms of the same number 
of independent quantities. This is because when we move one reaction from the first class 
to the second one, we lose one constraint from chemical equilibrium dynamics, but we 
gain one constraint from a new conservation law. While this is in general true, as we 
will see the detailed way in which this 'switching' of conditions is implemented has some 
subtleties, since in some cases a non trivial interplay between redundant conditions and 
new global symmetries that are anomalous, and thus do not correspond to conservation 
laws, is involved. 

This section is organized in the following way: first we list in Section |2.1| the con- 
straints that hold independently of assuming a regime in which particles and sparticles 
chemical potentials equilibrate (superequilibration (SE) regime) or do not equilibrate (non- 
superequilibration (NSE) regime). In Section we list the constraints that hold only in 
the SE regime, and in Section |2.3| the ones that hold in the temperature range relevant for 
super symmetric leptogenesis, when the temperature is sufficiently high that particles and 
sparticles chemical potentials do not equilibrate. 

2.1 General constraints 

The supersymmetric seesaw model is described by the superpotential of the super symmetric 
SM with the additional terms: 

W = ^MpgN;;N^ + Xp^N^ i^H^, (2.2) 

where p,q = 1,2,... label the heavy singlet states in order of increasing mass, and 
a = e,iJ,T labels the lepton flavour. In Eq. (|2.2[) £, Hu and A'''^, are respectively the 
chiral superfields for the lepton and the up-type Higgs SU (2) doublets and for the heavy 
SU{2) singlet neutrinos defined according to usual conventions in terms of their left-handed 
Weyl spinor components (for example the N'^ supermultiplet has scalar component N* and 
fermion component Nf^). Finally the SU{2) index contraction is defined as iaHu = epaiaH^ 
with ei2 = +1. 

Let us start at T ^ few TeV, that is well above the EW or SUSY breaking scale, but 
low enough that all the SM and SUSY reactions can be considered in thermal equilibrium. 
We first list in items 1., 2. and 3. the conditions that hold in the whole temperature range 
that we will consider {M]y <^ T ^ lO^^GeV). Conversely, some of the Yukawa coupling 

t Assigning B — L violating processes to class I means that they are in thermal equilibrium and no B — L 
asymmetry can be generated. Assigning them to class II means that B — L is conserved, and baryogenesis 
cannot occur. 
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conditions given in items 4. and 5. will have to be dropped as the temperature is increased 
and the corresponding reactions go out of equilibrium. For simplicity of notations, in 
the following we denote the chemical potentials with the same notation that labels the 
corresponding field: </> = fi^f,. 

1. At scales much higher than Myy, gauge fields have vanishing chemical potential W = 
B = g = [22]. This also implies that all the particles belonging to the same SU{2) 
or SU{3) multiplets have the same chemical potential. For example </>(/3 = +^) = 
0(/3 = — ^) for a field (f> that is a doublet of weak isospin /, and similarly for color. 

2. Fast gluino, wino and bino interactions relate the difference between the chemical 
potentials of the members of a given supermultiplet and the corresponding gaugino 
chemical potential. Furthermore Q + gR ^ Q, Q + Wr Q, £ + Wr i, i + 
Br — )• i, where i, Q {£, Q) denote the (s)lepton and (s)quarks left-handed doublets, 
and Wr, Br and gR are respectively right-handed winos, binos and gluinos, enforce 
chemical equilibrium conditions which imply that all gauginos have the same chemical 
potential: 

-g = Q-Q = -W = £-e = -B, (2.3) 

where we have introduced W, B and g to denote the chemical potential of the left- 
handed gauginos. It follows that the chemical potentials of the SM particles are 
related to the chemical potential of their respective superpartners as 

Q,i=Q,£+~g (2.4) 
Hu,d = Hu,d + g (2.5) 
u,d,e = u,d,e — g. (2.6) 

The last relation, in which u,d,e = UR,dR,eR denote the R-handed SU{2) singlets, 
follows e.g. from 'uP^ = u^^ + g for the corresponding L-handed fields, together with 



= —ur and from the analogous relation for the SU{2) singlet squarks. Eqs.(2.4)- 
( |2.6D allow us to express all the chemical equilibrium relations in terms of the 18 
chemical potentials of the fermions (SM quarks and leptons, higgsinos and gauginos). 

Before EW symmetry breaking hypercharge is an exactly conserved quantity. There- 
fore for the total hypercharge of the Universe we have 

3^tot = X Anb + X An/ y/ = const, (2.7) 
b f 

where ybj denotes the hypercharge of the 6-bosons or /-fermions. Given that lepto- 
genesis aims to explain the origin of the matter-antimatter asymmetry, it is reason- 
able to assume as initial condition that at sufficiently high temperatures all particle- 
antiparticle density asymmetries vanish An^ = which implies J'tot = 0. By 
using the hypercharge conservation condition can be rewritten as: 



Mfe 56 + ^ X] dfVf = ^ (Qi + 2'Ui - c?j) - X + + Hu- Hd = 0, 



a 



(2.8) 



-6- 



where Qi denote the three quark doublets, Ui = u, c, t and di = d, s, b. In writing 
down Eq.(p.g|) we have used Eqs.( |2.4| )-( |2.6| ) to express the chemical potentials of the 



scalars in terms of those of their fermion partners. Given that Eqs.(2.4)-(2.6) hold in 



the whole temperature range that we will analyze, the same is true also for Eq.(p.^). 

Note that the hypercharge neutrality condition involves only the 17 chemical potentials 
of the fermionic components of the matter superfields since the gaugino chemical potential, 
that is involved in the substitutions (2.4)-(2.6), eventually cancels out. This is expected, 



since all gauginos have vanishing hypercharge; however, in counting the number of chemical 
potentials g must be included, yielding a total number of 18. 

4. When the reactions mediated by the leptons Yukawa couplings are faster than the 
Universe expansion rate H ~ 1.66^T^/Mp (where Mp is the Planck mass and 
5* = 228.75 in the MSSM), the following chemical equilibrium conditions are enforced: 

ia - ea + Hd + g = 0, {a = e,fi,T). (2.9) 

For a = e the corresponding Yukawa condition holds only as long as 

T < 10^(1 + tan^ (3) GeV, (2.10) 

when Yukawa reactions between the first generation left-handed SU (2) lepton dou- 
blets ie and the right-handed singlets e are faster than the expansion [23]. Note also 
that, as is discussed in refs. [24,25], if the temperature is not too low lepton flavour 
equilibration induced by off-diagonal slepton soft masses will not occur. We assume 
that this is the case, and thus we take the three £a to be independent quantities. 

5. Reactions mediated by the quarks Yukawa couplings enforce the following six chemical 
equilibrium conditions: 

Qi - Ui + Hu + g = 0, {ui = u,c,t), (2.11) 
Qi-di + Hd+~9 = 0, {di = d, s, b) . (2.12) 

The up-quark Yukawa coupling maintains chemical equilibrium between the left and 
right handed up-type quarks up to T ~ 2 • 10^ GeV. Note that when the Yukawa 
reactions of at least two families of quarks are in equilibrium, the mass basis is fixed 
for all the quarks and squarks. Intergeneration mixing then implies that family- 
changing charged-current transitions are also in equilibrium: b^ — )• cl and ti ^ sl 
imply Q2 = Q3; sl — )• ul and cl — )• dL imply Qi = Q2. Thus, up to temperatures 
T ^ 10^^ GeV, that are of the order of the charm Yukawa coupling equilibration 
temperature, the three quark doublets have the same chemical potential: 

Q = Q3 = Q2 = Qi. (2.13) 

At higher temperatures, when only the third family is in equilibrium, we will have 
instead Q = Q?j = Q2 ^ Qi- Above T ~ 10^^ when (for moderate values of tan /3) 
also the r and 6-quark SU{2) singlets decouple from their Yukawa reactions, all 
intergeneration mixings are also negligible and Qz Q2 Q\- 



-7- 



2.2 Superequilibration regime 

At relatively low temperatures, additional conditions from reactions in chemical equilibrium 
hold. Since the constraints below apply only in the SE regime, we number them including 
this label. 

6sE- Equilibration of the particle-sparticle chemical potentials /z,^ = pL^ (generally referred 
as superequilibration [26]) is ensured when reactions like U — >• U are faster than the 
Universe expansion rate. These reactions are induced by gaugino interactions, but 
since they require a gaugino chirality flip they turn out to be proportional to its mass 
rrig, and can be neglected in the limit rrig — t- 0. 

Furthermore, since the ^ parameter of the HuHd superpotential term is expected to 
be of the order of the soft gaugino masses, it is reasonable to consider in the same 
temperature range also the effect of the higgsino mixing term, which implies that the 
sum of the up- and down- higgsino chemical potentials vanishes. The rates of the 
corresponding reactions, given approximately by Tg ~ rn^/T and ~ l^'^/T, are 
faster than the Universe expansion rate up to temperatures 

The corresponding chemical equilibrium relations enforce the conditions: 

9 = 0, (2.15) 
Hu + Hd = 0. (2.16) 

7sE- Up to temperatures given by (|1|) the MSSM has the same global anomalies than 
the SM, that are the EW SU{2)-U{1)b+l mixed anomaly and the QCD chiral 
anomaly. They generate the effective operators Oew = ^aiQQQ^a) and Oqcd = 
Ili{QQu'i^d'i^) . Above the EW phase transition reactions induced by these operators 
are in thermal equilibrium, and the corresponding conditions read: 

9Q + ^^a = (2.17) 
6Q -Y,iui + di) = 0, (2.18) 

i 

where we have used the same chemical potential for the three quark doublets (Eq. ( |2.13| ) ) , 
which is always appropriate in the SE regime below the hmit (|2T|). 



2.2.1 Flavour charges 

Eqs.(|]|) and (|T1|)-(|1|), together with the SE conditions (|]l|)-(|2l|), the two anomaly 
conditions (|2.17| )-( |2ri8|) and the hypercharge neutrality condition ( |2.8D , give 11 + 2 + 2-1-1 = 
16 constraints for the 18 chemical potentials. Note however that there is one redundant con- 
straint, that we take to be the QCD sphaleron condition, since by summing up Eqs.([2Tll) 
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and ( 2.12 ) and taking into account ( 2.13| ), ( 2.151) , and ( 2.16 ) we obtain precisely Eq.( 2.1g| ). 
Therefore, like in the SM, we have three independent chemical potentials, that could be 
taken to be the ones corresponding to the leptons doublets. Another choice, that is more 
useful in leptogenesis, is to define three linear combinations of the chemical potentials cor- 
responding to the SU{2) anomaly free flavour charges = B /3 — L^- The reason is that 
these charges, being anomaly free and perturbatively conserved by the low energy MSSM 
Lagrangian, evolve slowly because the corresponding symmetries are violated only by the 
heavy Majorana neutrino dynamics. Their evolution is thus determined by reactions be- 
longing to class III, and needs to be computed by means of three independent Boltzmann 
equations. We define the number densities of particles per degree of freedom, normalized 
to the entropy density s, as ^Afe,A/ = terms of these quantities the density of 

the Aq charges normalized to the entropy density can be written as: 



(2.19) 



The expression above is completely general and holds in all temperature regimes, including 
the NSE regime (see Section |2.3| ). Note that g in the equation above cancels for the quarks 
but not for the leptons, and thus in the NSE, in which the gaugino chemical potential 
does not vanish, when the Y/\^ charges are expressed just in terms of the number density 
asymmetries of the fermions, Y/^g also contributes. 

In Eq.(gl9D we have left in clear some numerical factors: the overall factor of 3 adds the 
contributions of scalars (that is twice that of fermions) , the factor of 2 in front of the Y^q. 
and Y^£^ accounts for the SU{2) gauge multiplicity, while the color factor compensates 
against the the quark baryon number B = 1/3. 

The density asymmetries of the doublet leptons and higgsinos, that weight the washout 
terms in the Boltzmann equations, can now be expressed in terms of the anomaly free 
charges by means of the A matrix and C vectors introduced respectively in ref. [8] and 
ref. [11] that are defined as: 



Yae 



Y, 



(2.20) 



Here and in the following we will give results for the A and C matrices for the fermion 
states. We recall that in the SE regime the density asymmetry of a scalar boson that is in 
chemical equilibrium with its fermionic partner is given simply by lAfe = ^^a/ with the 
factor of 2 from statistics. 

2.2.2 All Yukawa reactions in equilibrium 

Assuming moderate values of tan /3, at temperatures below the limit in Eq.( |2.10| ) standard 
leptogenesis cannot be successful. However, this range can correspond to a low temperature 
windows in which soft leptogenesis [27-29] can successfully proceed, and therefore it is 
worth giving the results for A and C. They are: 

/-221 16 16 



A' 



1 

9 X 237 



v 



16 
16 



-221 
16 




C 



-c 



Hd 



237 



:i, 1, 1). (2.21) 
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Note that since in this regime the chemical potentials for the scalars and leptons degrees 
of freedom of each chiral multiplet equilibrate, the analogous results for Ia^^ + ^^'^ 
be obtained by simply multiplying the A matrix in Eq.( |2.21 ) by a factor of 3. This gives 



the same A matrix obtained in the non-supersymmetric case in the same regime (see e.g. 
eq.(4.13) in ref. [11]). The C matrix (multiplied by the same factor of 3) differs from the 
non-supersymmetric result by a factor 1/2. This is because after substituting = —H^ 
(see Eq.( 2.16| )) all the chemical potential conditions are formally the same than in the SM 



with Hu identified with the chemical potential of the scalar Higgs, but since C expresses 
the result for number densities, in the SM a factor of 2 from boson statistics appears for 
the SM Higgs. This agrees with the analysis in ref. [30], and is a general result that holds 
for supersymmetry within the SE regime. 

2.2.3 Electron and up-quark Yukawa reactions out of equilibrium 

Raising the temperature above 10^(1 + tan^ /3) GeV the interactions mediated by the elec- 



tron Yukawa he are not able to maintain equilibrium, and one condition in Eq.(2.9) for 
a = e is lost. However, since in the effective theory at this temperature one can set /ig — )• 0, 
one global symmetry is gained. This corresponds in the fermion sector to chiral symmetry 
for the R-handed electron, that in the present case translates into a symmetry under phase 
rotations of the e chiral multiplet that holds in the limit of unbroken supersymmetry. Con- 
servation of the corresponding charge ensures that Aug + Ang = SAUf, is constant, and 
since leptogenesis aims to explain dynamically the generation of a lepton asymmetries we 
set this constant to zero, so that the R-handed electron chemical potential is e = 0. In this 
way the chemical equilibrium condition that is lost is replaced by a new condition implied 
by the conservation of a global charge, and three independent chemical potentials (or al- 
ternatively the three non- anomalous charges ( ^.19 )) are still sufficient to describe all the 



density asymmetries of the thermodynamic system. At temperatures above T ~ 2-10® GeV 
interactions mediated by the up-quark Yukawa coupling drop out of equilibrium. In this 
case however, by setting /i„ — )• no new symmetry is obtained, since chiral symmetry for 
the R-handed quarks is anomalous and the corresponding charge is not conserved by fast 
QCD sphaleron interactions. However, after dropping the first condition in Eq.( |2.11 ) for 
Ui = u, the QCD sphaleron condition Eq.( |2.18| ) ceases to be a redundant constraint, with 
the result that also in this case no new chemical potentials are needed to determine all the 
particle density asymmetries. The corresponding results, that can be again relevant for 
soft leptogenesis, read: 



C^" = -C^d = (37, 52, 52) .(2.22) 

2886 ^ ' ' ^ ^ ^ 



2.2.4 First generation Yukawa reactions out of equilibrium (SE regime) 

Let us now consider what happens at temperatures T > 4 • 10^(1 + tan^ /3) GeV, when 
also the d-quark Yukawa coupling can be set to zero (in order to remain within the SE 
regime we assume tan /3 ~ 1). In this case the equilibrium dynamics is symmetric under 
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the exchange u ■v^ d (both chemical potentials enter only the QCD sphaleron condition 
Eg. ( p. 18 ) with equal weights) and so must be any physical solution of the set of constraints. 



Thus, the first condition in Eq.( |2.12 ) can be replaced by the condition d = u, and again 



three independent quantities suffice to determine all the particle density asymmetries. The 
corresponding result is: 

/-906 120 120 

75 -688 28 I , C^" = -C^'' = — ^ (37, 52, 52) , (2.23) 



A 



£ 



1 



3 X 2148 



V 



75 28 -688 



2148 



that agree with what is obtained in non-super symmetric leptogenesis (see eq. (4.12) of 
ref. [11]) after the factor 1/2 for the higgsinos discussed below Eq. ( p. 211 ) is accounted for. 
In the numerical analysis in Section ^ we will take this case as a benchmark to confront 
the results obtained in the SE regime with the corresponding results in the NSE regime. 

2.3 Non-superequilibration Regime 

At temperatures above the limit given in Eq.( p.l4 ) the Universe expansion is fast enough 
that reactions induced by mg and fi do not occur. Setting to zero in the high temperature 
effective theory these two parameters has the following consequences: 

(i) Condition ( p.l5| ) has to be dropped, and gauginos acquire a non-vanishing chemi- 
cal potential g ^ (corresponding to the difference between the number of L and R 
helicity states). The chemical potentials of the members of the same matter supermul- 



tiplets are no more equal (non-superequilibration) but related as in Eqs.(2.4)-(2.6) 



ii) Condition ( p. 16 ) also has to be dropped, and the chemical potentials of the up- and 



down-type Higgs and higgsinos do not necessarily sum up to zero. 

(iii) The MSSM gains two new global symmetries: nig — )• yields a global i?-symmetry, 
while ;U — )• corresponds to a global symmetry of the Peccei-Quinn (PQ) type.-f 

2.3.1 Anomalous and non-anomalous symmetries 

The charges of the various states under the R and PQ symmetries together with the values 
of the other two global symmetries B and L are given in Table |. Like L, also R and PQ 
are not symmetries of the seesaw superpotential terms MN^N^ -\- XN^lHu- In Table |l] we 
have fixed the charges of the heavy A^'^ supermultiplets in such a way that the mass term 
MN'^N'^ is invariant under all symmetries and in particular has R charge equal to 2.^ In 

'"We assume that, similarly to rrig, in this regime also the other soft supersymmetry breaking terms 
can be neglected, and in particular the B-term for the heavy sneutrinos. This is always the case in the 
high temperature regime for successful leptogenesis. However, at the lower temperatures required in soft 
leptogenesis this is not necessarily true. Then the ii-charge is not perturbatively conserved, as is clear 
from the last entry in Table ^. This can have quite interesting consequences for soft leptogenesis that are 
analyzed in a companion paper [31]. 

^ Under i?-symmetry the superspace Grassmann parameter transform as ^ ^ e^°'9 . Invariance oi J d9d — 
1 then requires R{d6) = — 1. Then the chiral superspace integral of the superpotential J d6^ W is invariant 
if R{W) = 2. By expanding a chiral supermultiplet in powers of 9 it is also clear that its R charge equals 
the charge of the scalar boson R{b) — R{f) + 1. 
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9 


Q 




d" 


e 




Hd 


Hu 


N" 


B 





1 

3 


1 

3 


1 

3 

















L 














1 


-1 











PQ 








-2 


1 


-1 


2 


-1 


2 





R 1 

b 


1 


-1 


-3 


1 


-1 


1 


-1 


3 





2 





-2 


2 





2 





4 


1 



Table 1: B, L, PQ and R charges for the particle supermultiplets that are labeled in the top 
row by their L- handed fermion component. Note that we use chemical potentials for the R- handed 
SU{2) singlet fields u, d, e that have opposite charges with respect to the ones for u'^, d'^, e*^ given 
in the table. The i?-charges for bosons are determined by R{b) — R{f) + 1. 



contrast, the Yukawa term N^LH^ violates L, PQ and R. All the four global symmetries 
B, L, PQ and R have mixed gauge anomalies with SU{2), and R and PQ have also mixed 
gauge anomalies with SU{3). Two linear combinations R2 and R^ of R and PQ, having 
respectively only SU{2) and SU{3) mixed anomalies have been identified in ref. [21]. They 
are: ^ 

R2 = R-2PQ (2.24) 
R3 = R-3PQ. (2.25) 

The values of i?2,3 for the different states are given in Table |2[ The authors of ref. [21] have 
also constructed the effective multi-fermions operators generated by the mixed anomalies: 

Oew = n„ {QQQta) HuHd , (2.26) 
Oqcd = n, {QQu'd% t . (2.27) 

Given that three global symmetries B, L and R2 have mixed SU{2) anomalies (but 
are free of SU{3) anomalies) it is clear that we can construct 2 anomaly free combina- 
tions, the first one of which we chose to be, as in the SM, B — L. Given that only one 
term, that is IHu, violates perturbatively the global symmetries, a second anomaly free 
combination IZ can be chosen in such a way that it is respected also by the corresponding 
interactions, and thus it is an exact symmetry of the MSSM+seesaw in the NSE regime. 
The corresponding charge reads 

n=^B-L + R2, (2.28) 

and is exactly conserved. In the SU (3) sector, besides the chiral anomaly we now have also 
i?3 mixed anomalies. Thus also in this case anomaly free combinations can be constructed, 

^For definiteness we restrict ourselves to the case of three generations Ng = 3 and one pair of Higgs 
doublets Nh = 1. We also normalize _R2,3 in such a way that R2,3{b) = R2,-i{.f) + 1- 
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and in particular we can define one combination for each quark superfield. Assigning to 
the L-handed supermultiplets chiral charge x = — 1 these combinations have the form: 

Xq, + 1^,^ Rs (2.29) 

where, for example, = n^c^ = 1/3 and kq^ = 2/3. Note that since R^ is perturbatively 
conserved by the complete MSSM+seesaw Lagrangian, when the Yukawa coupling of one 
quark is set to zero the corresponding charge Eq. ( 2.291 ) will be exactly conserved. 



2.3.2 Constraints in the non-superequilibration regime 

As explained in the introduction to this section, conditions corresponding to conservation 
laws constrain particle number density-asymmetries, while conditions from reactions that 
are in chemical equilibrium constrain particle chemical potentials. Eventually, to get a 
closed form solution to the set of constraining equations one needs to use a single set 
of variables. In the following, we will then express the number density- asymmetries of 
relativistic particles in terms of their chemical potentials by means of eq. (|2.1| ) that holds 
in the ultra-relativistic limit and at first order in ^j {,/T. 

In the NSE regime, the conditions listed in items 65^; and 7se of the previous section 
have to be dropped, but new conditions arise. 

^NSE- The conservation law for the TZ charge yields the following global neutrality condition: 



r2 



^^/(7/% + 2 I - An^^ =0, (2.30) 



6 

with the term in parenthesis given by: 
fJ-f 9f + 2 ^ fJ-b Qb Tib = 

f b 

2 (^(2Qi-5tx, + 4(i,) + 2^(^„ + e„) + 5^d-^„ + 3l5 ) . (2.31) 



The last terms in both lines of Eq. ( |2.30| ) correspond to the contribution to TZ- 
neutrality from the lightest sneutrino asymmetry An^^ = ~ '^tv* '^ith charge 
TZj^_^ = —TZ]\fc = —1. All the other heavier neutrinos are assumed to have already 
decayed at the time A^i leptogenesis takes place, and do not contribute. Note that 
since in general Ni is not in chemical equilibrium, no chemical potential can be asso- 
ciated to it, and hence this constraint needs to be formulated in terms of its number 
density asymmetry that has to be evaluated by solving a Boltzmann equation for 
YAj^ = Yff-^ — ^7v* (^^^ Section 3). Note also that since the final lepton asymmetry 
in leptogenesis is generally determined by the dynamics at T < M, where the sneu- 
trino abundance gets exponentially suppressed, we can expect that neglecting An^ 



in Eq. ( 2.30 ) could still give a good approximation. This is indeed confirmed by the 



numerical analysis presented in Section |^. 
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'^NSE- The operators in Eqs. ( 2^61 ) - (|2 .271 ) induce transitions that in the NSE regime are in 
chemical equilibrium. This enforces the generalized EW and QCD sphaleron equilib- 
rium conditions [21]: 



i a 



(2.32) 
(2.33) 



that replace (|2?T7|) and (pl|) 



8nse- The chiral-i?3 charges in Eq.( 2.29 ) are anomaly free, but clearly they are not con- 
served by the quarks Yukawa interactions. However, when a quark supermultiplet 
decouples from its Yukawa interactions an exact conservation law arises. (Note that 
> implies u and d decoupling, but Qi decoupling is ensured only if also 
0.) The conservation laws corresponding to these symmetries read: 



hu,d 



— [3gR -I- 6(gR -g)] + - R3 tot = 

^ 2 [3Ql + 6{Ql + ~9)]-l i?3tot = 
3 



(2.34) 
(2.35) 



and hold for qji = Ui, di and Qi = Qi in the regimes when the appropriate Yukawa 
reactions are negligible. Note the factor of 2 for the Ql chiral charge in front of the 
first square bracket in Eq. ( p.35| ) that is due to SU (2) gauge multiplicity. In terms 
of chemical potentials and of the sneutrino number density asymmetry, the total 
charge in Eqs. ( |2.34| )-( |2l3^ ) reads: 



R. 



'3 tot 




(2.36) 



where R^ 



-1, and the quantity in parenthesis is given by: 
^lif gf R3f + 2^fib9b R3 



I3b 



82g - 3^ (2Qi + llu, - 4di) + (16^« + 13e«) + 16 Hd - UHu- (2.37) 



As regards the leptons, since they do not couple to the QCD anomaly, by setting 
/le — >• a symmetry under chiral supermultiplet rotations is directly gained for the 
R-handed leptons implying Arie + Ang = and giving the condition: 




(2.38) 



No analogous condition arises for the lepton doublets relevant for leptogenesis, since 
by assumption they remain coupled via Yukawa couplings to the heavy A^'s. 
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Q 




(F 








Hu 




b 


1 


-1 


1 


-1 


1 


-3 


1 


-1 





2 





2 





2 


-2 


2 





1 


b 


1 


-1 


3 


-2 


2 


-5 


2 


-3 





2 





4 


-1 


3 


-4 


3 


-2 


1 


f 
b 


1 


4 

9 


4 
9 


14 
9 





-2 


1 


-1 





2 


5 
9 


13 
9 


5 
9 


1 


-1 


2 





1 



Table 2: Charges for the fermionic and bosonic components of the SUSY multiplets under the 
i?-symmetries defined in Eqs.(2.24), (2.25) and ( 2.28| ). Supermuhiplets are labeled in the top row 
by their L-handed fermion component. We use chemical potentials for the R-handed SU{2) singlet 
fields u, d, e that have opposite charges with respect to the ones for u'^, , e'^ given in the table. 



Let us now evaluate what is the number of independent chemical potentials in the NSE 
regime. Since we already know that this number should not change when some constraints 
from chemical equilibrium are dropped and replaced by conservation conditions, let us 
start by assuming tan j3 large enough that the Yukawa couplings he and hd enforce the 



equilibrium conditions (2.9) and (2.12), and also keep for the time being the condition 



( |2.11 ) for the up-quark. We have 9 conditions from the Yukawa couplings plus 2 from 



quark doublets equilibration ( p. 131 ), the generalized EW sphaleron condition (2.32) and 



the QCD sphaleron condition ( |2.33| ) that now must be counted since it is independent 
from other constraints, the global hypercharge ( |2.8| ) and IZ ( |2.31[ ) neutrality conditions, for 
a total of 15 constraints for the 18 chemical potentials or corresponding particle density 
asymmetries considered in the previous section. However, in the NSE regime we have that 
the set of constraining conditions involve one additional independent quantity that is the 
sneutrinos density asymmetry ^a^, and therefore the total number of relevant quantities is 
19. We thus conclude that in the NSE regime the three flavour charges Ya^ iii Eq.( 2.19|) do 



not suffice to describe the density asymmetries of all MSSM particles and superparticles, 
and that Ya^ is also required. 

To proceed, we can now set /i„ — )■ to match the actual number of equilibrium con- 
strains in the NSE regime. We lose one Yukawa condition, but we gain the corresponding 
chiral-i?3 conservation condition ( 2.34]) and, as expected, the overall number of independent 



chemical potentials or number density asymmetries remains the same. 

In the following we give some results corresponding to temperature ranges relevant for 
successful supersymmetric leptogenesis, that always occurs within the NSE regime. Given 
that the Yukawa conditions for the leptons and for the down-type quarks depend on the 
value of tan/3, the temperature ranges we refer to is only indicative (and correspond to 
moderate values of tan/3), but we will always specify in a clear way which conditions are 
used to obtain the corresponding results. In the non-superequilibration regime there are 
different flavour mixing matrices for the scalar and fermion components of the leptons 
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and Higgs supermultiplets. To express more concisely all the results, it is convenient to 
introduce a new C vector to describe the gaugino number density asymmetry per degree 
of freedom in terms of the relevant charges: 

Y^g = Cl , with A„ = (A„, Aj^) . (2.39) 

2.3.3 First generation Yukawa reactions out of equilibrium (NSE regime). 

In the temperature range between 10® and 10^^ GeV, and for moderate values of tan /3, all 
the first generation Yukawa couplings can be set to zero. Using for u, d conditions ( 2.34| ) 
and for e condition ( |2.3§| ) as are implied by /le — ^ we obtain: 

/ -198117 33987 33987 -8253 ' 
26634 -147571 14761 -8055 
26634 14761 -147571 -8055 



1 

9 X 162332 



= — (163, 165, 165, -255) , 

162332 ^ ' ' ' ^ ' 

C""- = (3918, 4713, 4713, 95) , 

C""' = 3 ^ ^g^33^ (5413, 9712, 9712, -252) , (2.40) 



where the rows correspond to (Ya^, Ya^, Ya^, ^a^^)- For completeness, in Eq. ( |2.4C ) we 
have also given the results for C^'' even if only the up-type Higgs density asymmetry is 
relevant for the leptogenesis processes. Note that neglecting the contribution of An^^ to 
the global charges T^tot in Eq. (|2.3C| ) and -Rstot in Eq. ( ^.36 ) corresponds precisely to set 



to zero the fourth column in all the previous matrices. Then, analogously with the SE and 
SM cases, within this '3-columns approximation' all particle density asymmetries can be 
expressed just in terms of the three Ya^ charge densities. 

In the numerical analysis carried out in the next section we will confront the results 
obtained with the full matrices ( p. 40 ) with the results obtained for the analogous case 



discussed in Section 2.2.4 of first generation Yukawas out of equilibrium, but within the 



SE regime, that yielded the matrices in ( |2.23| ). 

2.3.4 Second generation Yukawa reactions out of equilibrium: two flavour regime. 

In the temperature range between 10^^ and 10^^ Q^y 

we set hfj^ — t- 0. This condition defines 
the two flavour regime for leptogenesis [8,10,11]. In this regime we have also hg, he ^ that 
result in conserved Xql + i^gl^s charges (g = s, c) yielding the corresponding conditions 
( |2.34 ). Since the first generation quark doublet does not mix with the other generations. 



Qi is not equal to Q = Q2 = Qs, but conservation of XQi + i^Qi^s then provides the 
required new condition ( p. 35 ). Conservation laws stemming from h(,^^^u,d,s,c ~^ then add 



up to a total of 7 constraints that replace the 6 Yukawa equilibrium conditions plus the 
condition Qi = Q. 

In the two flavour regime only the lepton doublet ir is identified, while the ig and 
doublets, both of which have no Yukawa interactions, are not distinguished. In general. 
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one combination i^^ of these two doublets remains coupled to lightest heavy singlet A^i, 
while the orthogonal combination decouples, so that A^i decays do not generate directly a 
lepton asymmetry in this particular flavour direction. The corresponding non-anomalous 
B/3 — La charge is thus exactly conserved, and we can set its value to zero. This condition, 
together with the third generation Yukawa conditions, the two sphaleron conditions, and 
y and TZ global neutrality, give 8 constraints for the remaining 10 chemical potentials 
Q, b, t, £a, T, Hu4, g- Therefore, in this case the two flavour charge densities Iae^ and Ya^ 
together with Y/\^ suffice to express the density asymmetries for all the other particles. In 
the basis {Y^^^, Ya^, ^a^v) we obtain: 



A' 



1 / -460047 88166 -17229 \ 



6 X 580163 \ 80783 -349740 -14094 / 



C~3 = (-11056, -11786, 20307) , 

2 X 580163 ^ ^ 

= 27^5163 ™''''' 

C^<* = (8939, 77012, -12483) . (2.41) 

6 X 580163 ^ ' ' ^ ^ ^ 



Note that since EW sphaleron flavour mixing involves all the three lepton doublets, a 
non- vanishing lepton asymmetry (equal to B/3) is induced also in the decoupled flavour 
direction orthogonal to £f>^. However, since the corresponding doublets do not participate 
in the leptogenesis dynamics, this asymmetry is irrelevant for baryogenesis, and thus giving 
just the two entries corresponding to ir and ie/i in the matrices A and C above is sufficient. 



2.3.5 Above the EW sphaleron equilibration temperature: the B = regime. 

EW sphaleron processes take place at a rate per unit volume T/V oc T^a^ log(l/avK) 
[32-35], and are expected to be in equilibrium up to temperatures of about ~ 10^^ GeV. 
Even for moderate values of tan/3, the b and r Yukawa interactions are still in equilibrium 
at this temperature. We will then consider here the temperature regime 10^^ — 10^^ GeV 
in which super symmetric leptogenesis still occurs in the two flavour regime, but with EW 
sphalerons switched off. The related chemical equilibrium condition ( |2.32| ) is then replaced 
by baryon number conservation. We set i? = as initial condition since a dynamical 
generation of the baryon asymmetry is precisely the goal of leptogenesis. This also implies 
that Ya^ = —YaLci and therefore the total lepton asymmetry in the direction orthogonal 
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to defi vanishes exactly. For this regime we obtain: 

1 / -72807 1480 -5676 \ 
6 X 73327 i -77 -50984 -4236 j ' 



A' 



= (-130, -370, 1419) , (2.42) 

= ^ ^~33^y (9187, 9226, 686) , (2.43) 

= a ^ 7^-^97 (1^2^' 9998, -1482) . (2.44) 
6 X 73327 



2.3.6 Above the temperature for r Yukawa equilibration: the one flavour regime. 

For moderate values of tan /3, at temperatures above T ~ 10^^ GeV we can set hr ^ 0. 
This condition defines the one flavour regime, in which only the dynamics of the lepton 
doublet ii to which A^i couples is important. This allows us to set to zero the two lepton 
density asymmetries in the direction orthogonal to £i : YaLi = Yal' = l^ft 
with just one non- vanishing lepton flavour asymmetry that is Y^Li- On the other hand, 
approximate b — t Yukawa unification suggests that we should also set hh — )■ 0, so that 



only the ht Yukawa condition (2.11) remains. Since quark mixing between the second and 
third generation is, for small tan/3, of the same order than the ratio h^/ht ~ 0(10^^), 
we assume that Q2 is also decoupled from Yukawa interactions and thus we replace the 
condition Q2 = Q3 by the conservation of the corresponding charge XQ2 + I'^Qi^s^ ^^^^ 
yields the constraint ( 2.35| ) . 



Note that even if in this regime all lepton flavour effects can be neglected [8, 10, 11], 
not all spectator effects can [7]. For example, QCD sphaleron equilibration is maintained 
up to higher temperatures than the corresponding electroweak processes because they have 
larger rates: ^qcd/^ew — 12(as/aiy)^, and are likely to be still in equilibrium at T ~ 
10^^ GeV [35-37]. Since Y^.^ remains coupled to the lepton flavour dynamics because of 
its contribution to T^tot and -Ratot we still have matrices with two columns (— Yal^i ^a^^^) 
even in the one flavour regime. After imposing all the relevant conservation conditions, 
including the vanishing of lepton number in two flavour directions, we obtain: 

A' = -^(49, 4). CS = ^(-1,24), 

Before concluding this section, it is worth reminding that the conversion of the B — L asym- 
metry generated by the leptogenesis dynamics into a baryon asymmetry should eventually 
be computed at a temperature 0(100 GeV), that is right before EW sphaleron transitions 
are switched off. At this temperature presumably all the supersymmetric particles have 
already decayed, and the particle content of the thermal bath is then the same than in the 
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SM with two Higgs doublets. Depending if this temperature is higher or lower than the 
temperature Tewpt of the EW phase transition, the conversion factors are [38]: 



Yab = lA._. X ^ - T > Wt (2.46) 
3. Numerical Results 

As we have seen in the previous section, supersymmetric leptogenesis is characterized by 
important qualitative differences with respect to SM leptogenesis. These differences arise 
because in the NSE regime, which is the relevant one to ensure that a sufficient amount 
of baryon asymmetry can be generated, supersymmetric leptogenesis cannot be treated by 
simply accounting for the new degrees of freedom in the leptogenesis dynamics. In this 
section we analyze the quantitative relevance of the new effects we have been discussing. 
They are all related with washout effects that in the NSE regime are controlled by the 3x4 
matrix Wnse = ^nse + ^nse + ^nse^ while in the SE regime the corresponding matrix 
yVsE = + Cg^ is 3 X 3. We can distinguish three types of effects: 

(i) The overall strength of the washout processes, that is mainly controlled by the diag- 
onal entries of the 3x3 submatrix W^^^, has slightly larger weights. 

(ii) Lepton flavour mixing effects, that are controlled by the off diagonal entries in Wf^g^ 
are sizeably larger in the NSE regime. 

(iii) New contributions to the washout arise from mixing with ^A^^;,) given by the entries 
in the fourth column W^'g^- 

One can get an intuitive handle about the quantitative impact of these effects on the final 
result for the baryon asymmetry, by comparing the overall washout coefficients evaluated 
in the SE regime for the case of the first generation Yukawa reactions out of equilibrium 
given in Eq. ( p.23| ), to the corresponding coefficients in the NSE regime given in Eq. ( 2.4C| ). 
We find, approximately independently of the particular entry: 

i^Nsl - ^se]^^ ^ « -0.01. (3.1) 

For the (larger) diagonal entries this difference remains somewhat below 10%; for the 
(smaller) off-diagonal entries in W^^^^; this represents an 0(1) difference, while the W^^^; 
entries are specific of the NSE regime. We have found that in the regimes in which the 
washouts are rather strong, and depending on the specific fiavour configuration, the first 
two effects can produce 0(1) changes in the final value of the baryon asymmetry. The last 
correction, although potentially of a similar size, gives some effects on the evolution of the 
flavour density asymmetries only for z ^ 1, that is when the A^'s approach an equilibrium 
distribution. However, at z > 1, Y/\^ gets exponentially suppressed implying that basically 
no effects are left in the final result. 

The results of our numerical analysis are summarized by the plots in figures I and 0, 
that are obtained by integrating numerically the complete set of Boltzmann equations given 
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in the Appendix. They include decays, inverse decays and scatterings with top-quarks, and 
hold under the assumption that the A^'s are hierarchical and that the lepton asymmetries 
only result from the decays of the lightest heavy singlet states N ~ Ni. However, in order 
to illustrate how the new effects described above modify the structure of the equations, here 
we write much simpler expressions in which only decays and inverse decays are included: 



Yn 



Y, 



. N 



Ya^ 



N 




7iv 



Yn 



5^ 



ye? - 1 I 77V + 



N 



Y, 



N 



2 7^ 



+ 7| + ^7^ E Ka + C^-'+C! 



Ya^ 



(3.2) 
(3.3) 
(3.4) 



(3.5) 



In the equations above the time derivative is defined as 1" = zsH dY/ dz with s the entropy 
density, z = M/T and H = H{M) the expansion parameter at T = M. In Eq. (3.3) we 



have introduced the overall sneutrino abundance 



Y^ + Yj^,, while Fa^ 



Y, 



N* 



in Eq. (3^) is the sneutrino density asymmetry that was already introduced in Section 
2.3.2| . In the washout terms we have normalized the charge densities = (^Aa) ^A^) 
to the equilibrium density of a fermion with one degree of freedom Y^, and we refer to the 
Appendix for detailed definitions of the reaction densities and associated quantities. In 
Eqs. (^.2|)- (|3.5| ) we have also neglected for simplicity all finite temperature effects. Taking 
these effects into account would imply for example that the CP asymmetry for N decays 
into fermions is different from the one for decays into scalars, while we describe both CP 



asymmetries with ea- A few remarks regarding Eq. (3.4) are in order. In the SE regime 
^ = and thus it would seem that the sneutrino density asymmetry Ya^ vanishes. However, 
this only happens for decays and inverse decays, and it is no more true when additional 
terms related to scattering processes, that are represented in the equation by the dots, 
are also included (see ref. [18] and the Appendix). Therefore, also in the SE regime Y^^ 
and Yfjt in general differ. However, in this case recasting their equations in terms of two 
equations for Y^^ and Ya^ is just a convenient parametrization. On the contrary, in the 
NSE regime this is mandatory, because the sneutrinos carry a globally conserved i2-charge 
and Ya^ is required to formulate properly the corresponding conservation law. As we have 
seen, this eventually results in Ya^ contributing to the expressions of the lepton flavour 
density asymmetries in terms of slowly varying quantities. 

The results for the final baryon asymmetry (normalized to the total CP asymmetry 
e = ea) are depicted in Figure |. We have fixed the value of the N mass to M = 
10^" GeV that is, we work in the three flavour regime in which all the Yukawa couplings 
of the fermions of the first generation can be set to zero and the results of Section |2.3.3 
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Figure I: The baryon asymmetry normalized to tlie total CP asymmetry IYab/^I computed in the 



three flavour NSE regime with the A and C matrices in (2.40) (solid red lines) confronted with 



what is obtained in the SE regime ( ^.23 ) (dashed blue lines). The two upper panels (a) and (b) 



correspond to strong washouts with rhi =0.1 eV, the two lower panels (c) and (d) to weak washouts 
with nil = 5 X 10^^ eV. In the two panels (a) and (c) on the left a flavour equipartition configuration 
has been assumed, with flavoured parameters ec/e ^ Pa = 1/3. In the two panels (b) and (d) on 
the right a strongly misaligned flavour configuration has been used, with — 0, — —2er — e 
and Pe = 1/100, = 9/100 P^ = 90/100. 



hold. We confront our results with the SE case discussed in Section 2.2.4 in which the SE 
conditions 5 = and Ya^.aHu = "^^aI AH ^o^d. As it is customary we parametrize the 
overall strength of the washouts in terms of the effective mass [39,40] 

a a a 

where Vu is the vacuum expectation value of the up-type Higgs doublet, Vu = v sin (3 (f =174 
GeV).ll For the single lepton flavours, the corresponding parameters are defined in terms 
of the tree level N and N decays branching ratios Pa = |AiqP/(AA^)ii according to: 

rhia = Parhi, ^Pq, = 1. (3.7) 



" Once the set of Yukawa conditions used to derive the A and C matrices is established, the dominant 
dependence on tan j3 in the numerical results arises via Vu in Eq. (|3.6|) and therefore it is very mild. 
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In Figure | the solid red lines represent the results obtained with the NSE A and C 
matrices in Eq.( p.4C) ), while the dashed blue lines give the corresponding results obtained 
assuming SE and using A and C in Eq.( 2.23D . In the two upper panels (a) and (b) we 
give the results obtained assuming a particularly strong washout regime corresponding to 
nil = 0.1 eV, while the two lower panels (c) and (d) correspond to a weak washout regime 
with nil = 5 X 10~^ eV. In the attempt to disentangle the differences between NSE and SE 
due to the changes in the overall washout strength [effects of type (i)] from those due to 
the changes in the flavour mixing pattern [effects of type (ii)] we have assumed for the two 
panels (a) and (c) on the left a flavour equipartition configuration [^ = i; Pq = i] while 
for the two panels (b) and (d) on the right we have assumed a strongly misaligned flavour 
configuration [^ = (0, 2, —1); Pa = (jgo' TOo' Mj)]' ^^^^ ^^^^ differences 

in the fiavour mixing patterns produce irrelevant effects, and the changes that can be 
seen are essentially due to the different fiavour diagonal washouts. In the second case 
the differences in the patterns of fiavour mixing also contribute. As it is apparent from 
figure |, the numerical differences between the NSE and SE cases remain typically at the 
0(1) level. As regards the effects of type (iii) that are related to ^Aj^, we have verified that 
in all the plots in figure | any visible difference would appear only around z ^ 1, leaving 
the final results unchanged. However, for completeness, we present in figure || the detailed 
evolution of the single fiavour charge densities Ya„ and also of Ya^^ corresponding to the 
flavour /washout configuration of figure |(b). As in the previous figure, the NSE results are 
depicted with solid lines [from up to down and around z = 0.5: Ya^ (blue), Ya^ (magenta), 
Ya^ (black), Yab (red)] and the SE results are depicted with dashed lines. The thicker 
portion of each line corresponds to positive values of the corresponding asymmetry, while 
the thinner portion to negative values. A couple of interesting features can be identified: 
firstly, the largest differences occur for Ya^ a-^d this is because given that eg = 0, this 
asymmetry evolves solely because of flavour mixing effects; secondly, it can be seen how 
Ya^ gets strongly suppressed around the time when the fiavour asymmetries 'freeze out' 
(z ~ 5), and this explains the irrelevance of its effects on the final result. One can also 
notice that curiously in the NSE regime Ya^ changes sign twice, while only once in SE. 
This is, of course, due to the related mixing effects that are absent in the SE case. 



4. Discussion and Conclusions 



Motivated by all the recent advancements in leptogenesis studies, we have revisited the 
theory of supersymmetric leptogenesis in the attempt to put it on a more firm theoret- 
ical ground. By digging in some depth into the analysis of the various conditions that 
constrain the density asymmetries of the different particles in the thermal bath, we have 
found that important qualitative differences exist with respect to SM leptogenesis. We 
have first clarified the reasons why, even if the specific constraining conditions are differ- 
ent for the different temperature ranges in which leptogenesis can occur, the whole set 
of particle density asymmetries can be always expressed in terms of the same number of 
independent quantities. This fact is explained in the following way: whenever it happens 
that by raising the temperature a chemical equilibrium condition ceases to hold because 
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Figure II: Evolution of tlic normalized single density asymmetries |lA^/e| and |yAji;j/e| for the 
flavour/washout configuration of figure |(b). From up to down and around z — 0.5: Ya,, (blue), 
(magenta), Ya^^ (black) and Fa^ (red). Solid continuous lines give the results for the NSE 
regime, with the thicker portion of each line corresponding to positive values of the asymmetry, and 
the thinner portion to negative values. The SE results are depicted with dashed lines and with the 
same color/flavour correspondence. 



the particle reaction enforcing it goes out of equilibrium, then one Lagrangian parameter 
related to this out-of-equilibrium reaction can be set to zero. This generically results in 
a new global symmetry, and the associated conservation law then enforces a new condi- 
tion for the particle number densities that replaces the one from chemical equilibrium. 
Following this scheme, we have identified new symmetries that are specific of the super- 
symmetric leptogenesis case. New chemical potentials associated to the superpartners of 
the SM particles are also a generic feature of super symmetric leptogenesis, and a quite 
crucial example is provided by the gauginos. At the typical high temperatures relevant for 
leptogenesis, chirality hipping transitions for the gauginos are completely out of thermal 
equilibrium and basically do not occur. Gauginos thus develop a chemical potential that 
can be thought of as the difference between the number of L- and R-handed states. Three 
main consequences follow from this. The first one is that particle and sparticle chemical 
potentials are different during leptogenesis. This is because only at relatively low tempera- 
tures their chemical potentials are equilibrated by scatterings. However, the rates of these 
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scatterings vanish in the limit of vanishing gaugino masses, and are completely negligible 
at the high temperatures relevant for leptogenesis. This modifies the weights of the scalars 
and fermions-related washouts. The second one is that by setting to zero the gaugino 
mass, a new i?-symmetry arises. While this symmetry is anomalous, a suitable anomaly 
free combination can be constructed and this corresponds to an exactly conserved global 
charge. The related conservation constraint, that is specific to supersymmetry, modifies 
the A and C matrices that describe the mixing between lepton flavours induced by the 
EW sphalerons. The third one is that since sneutrinos carry a unit charge under the R 
symmetry, their number density asymmetry enters the corresponding conservation law. 
The result is that through the A and C matrices the lepton flavour asymmetries are now 
related to four independent quantities rather than the usual three flavour charges of the 
SM case, being the sneutrino density asymmetry the fourth one. 

While we believe that the results presented in this paper are quite interesting from the 
theoretical point of view, by comparing the baryon asymmetry yield of leptogenesis within 
our framework with the results obtained by neglecting all the new effects, we have concluded 
that quantitatively we are dealing with 0{1) corrections. This is hardly surprising since, 
quite in general, sizable changes in the baryon asymmetry generated through leptogenesis 
are related with the identification of new sources of CP violation or by spoiling some exact 
cancellation. While this is the case e.g. for fiavoured leptogenesis [8, 10, 11] that requires 
the inclusion of the fiavour CP asymmetries e^, or for soft leptogenesis [27-29] where new 
CP violating phases play a crucial role and thermal effects spoil an exact zero-temperature 
cancellation between the CP asymmetries for decays into scalars and fermions, in our case 
there are no changes in the amount of CP violation that drives leptogenesis nor cancellations 
that could be spoiled. 
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A. Boltzmann equations for supersymmetric leptogenesis 

The Boltzmann equations that we have used to derive the numerical results read: 
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Besides the decays and inverse decays included in Eqs. ( |3.2| )-( |3^ , these equations also 
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include scatterings with the top-quark, both in the washout and in the CP asymmetries. 
In the hmit in which neutrinos are sufficiently hierarchical, as we are assuming here, the 
CP asymmetries in scatterings with top quarks are the same than the CP asymmetries 
in decays [12, 13] and then can be easily included. We have not included gauge boson 
scatterings nor the corresponding CP asymmetries. We have neglected three body decays 
since their contribution is phase space suppressed, and also because we have found that in 
ref. [18], on which our equations are based, some diagrams related to these processes have 
been overlooked. AL = 2 scatterings mediated by off-shell singlet neutrinos are also left 
out since in the temperature range T ~ 10^*^ GeV in which our results are obtained they 
are completely irrelevant. We have also neglected all finite temperature effects except for 
the Higgs thermal mass that is kept to regulate the infrared divergences in scatterings with 
the top-quark. 

We define the equilibrium densities per degree of freedom normalized to the entropy 
density as: 

Y^" = = -y^" = Y'J = (A.5) 

with the MSSM number of effective degrees of freedom g^, = 228.75. The number density 
asymmetries are defined according to 1a£ — ~ ^I- Density asymmetries and reaction 
densities without an index a are understood to be summed over flavours: Ya^ = Yla ^A^„ 
and 7iv = 1^Q,7^'• In Eqs. (^^)-( |3^ ) the density asymmetries for scalars Y^j and 1ah„ 
can be expressed in terms of the corresponding asymmetries for fermions ^A^a and Y^jj 
according to: 

Veq ^eq "T ^eq j l^-UJ 

where Y-" = Y~^ = y/''. In the SE regime, Ya~ = and thus ^ = ^ and - 



Y - 



follow. The reaction densities entering Eqs. ( |A.1| )-(A.4) are: 
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The reduced cross sections for the processes hsted above can be found in ref. [18]. 
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